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Abstract 

Methods  are  described  for  ceJ-culatlng  the  energy  levels  for 
the  overall  rotation  and  internal  torsion  of  molecules  consisting 
of  a rigid  symmetrical  top  attached  to  a rigid  asymmetrical  frame - 
vork  in  such  a way  that  the  symmetry  axis  of  the  top  coincides  with 
a principal  axis  of  the  molecule.  Probable  exasqples  are  nitromethane 
and  CH^Fg.  Matrix  perturbation  methods  are  eniployed  to  obtain 
finite  rotational  secular  equations  veLlld  in  each  of  the  coses: 
lov  barrier,  high  barrier,  lov  asymmetry.  These  secular  equations 
are  modifications  of  the  ordinary  Wang  equation  for  the  rigid 
asymmetric  rotor  and  can  usiuiUy  be  solved  by  the  continued  fraction 
method.  The  symmetry  groves  applicable  to  this  problem  are  e^.so 
discussed. 
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Molecules  sucli  as  or  nl  tron.etlA.T-e'^-  (OE^SC^'  shew  overall  rotation - 

internal  torsion  of  the  methyl  group  relati.ve  to  the  groiQi,  &ni  vibration. 

1 

Microwave  spectra  due  to  these  mr*tions  h&ve  been  recently  reportea  for 
nitromethane  by  Tanneribaum^  Johns  on Myers,  and  lwln?z.  J.  Chem.  Phys. 

Sk9  (195^). 

For  the  purpose  of  treating  the  overall  rotation  and  ivitemal  torsion,  it  is 
convenient  to  think  of  the  BFg  or  SCg  grens*  is-s  a ri  Jiia  framework  to  which  is 
attached  a symmetrical  top  (CH^)  which  n&"  rotate  or  v. If  rate  about  an  axis 
colinear  with  a principal  axis  of  inertia  of  the  wl:ole  molecule.  Purtbennore, 
the  three  principal  moments  of  inertia  cf  the  whole  morec-fie  are  all  different 
and  are  uninfluenced  by  the  intemal.  orientstior.  cf  tho  ^tt  icbed  top.  The 
calc\xlation  of  the  rotational  ani  internal  torsicurl  e..cr®”  levels  of  this 

r 

class  of  molecules  is  the  subject  of  the  present  ptvVAV.“' 

2 

The  class  treated  here  is  in  some  wajB  simpler  .nan.  the. case  of  methc^l 
alcohol  which  has  been  thoroughly  studied  by  I)enni.3Cc;  .and,  ccwc<rkere  (see 
Ivash  and  Dennison,  J.  Chem.  Phys.  gl,  .l8C*+  (r.95^)  } but  the  approa.cb  differs 
somewhat. 

Model  and  Coordinates 

The  model  then  consists  of  two  connected  rigid^  bodies,  one  (the  top) 

^ For  a dlscTisslcn  of  vibrations  in  molecules  of  this  type,  see  B.  L. 
Crawford,  Jr.  and  E.  B.  Wilson,  Jr.,  J.  Chem.  Phys.  9.»  523  (19^1). 

having  two  equal  principal  moments  of  inertia  about  axes  perpendicular  to  a 
principal  axis  of  the  whole  moleciLle.  There  are  four  degrees  of  freedom, 
three  for  overall  rotation  and  one  for  rotation  of  the  top  about  its  unique 
axis.  See  Fig.  1. 


Lot  XfjfZ  "be  moving  cartesian  axes  rigidly  attached  to  the  framework  part 
of  the  molecule  and  coincident  wj.th  the  principal  axes  of  inertia  of  the  whole 
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molecule  (origin  is  at  center  of  mass  of  whole  moleciile).  The  z axis  will 
coincide  with  the  symmetry  axis  of  the  top.  The  three  Valerian  angles  9,^, 
and^C  of  x^y^z  relative  to  space-fixed  axes  will  describe  the  overall  orientation 
of  the  molecule,  while  the  angle  ct  will  give  the  relative  orientation  of  top 
and  framework. 

The  Hamiltonian 

k 

The  kinetic  energy  has  previously  been  written  as 
B.  L.  Crawford,  Jr.,  J.  Chem,  Phys.  8^,  275 

ar  - + ij?  * (i) 

in  which  is  the  moment  of  Inertia  of  the  top  about  its  symmetry  eocls, 

I are  the  principal  momenta  of  the  entire  rooi  >cule,  and  U) 

X jT  z X y 

sre  conqpcnenta  of  angular  velocity  of  the  framework  along  x,y,z.  It  has 
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been  uaual*^  to  eliminate  the  cross-product  terms  by  means  of  a trsms formation, 

^ iiee  K.  S.  Pltzer  and  W.  B.  Gwinn,  J.  Chem.  Phys.  428  (1942). 

but  this  will  not  be  done  in  this  treatment. 

To  obtain  the  Hamiltonian  form,  use  the  definitions  of  the  momenta: 

p * ^T/da,  « d’S/dCj^f  etc.  (2) 

This  leads  in  the  usual  way  to  the  form 

H - AP  ^ + BP^^  + CP  ^ - 2CpP,  + V(a) 

X y z z 


(5) 


k 


in  which  all  smgular  momenta  are  In  units  = h/27T*  and 


m 


V(a)  is  the  potential  energy  restricting  the  internal  rotation.  In  the 

model  \ised  here,  the  coefficients  are  all  constants.  Note  that  C Involves 

the  mojnent  of  Inertia  of  the  framework  part  alone  while  F contains  the  reduced 

moment  of  the  two  parts  of  the  molecule. 

I<he  quantities  F , P and  P were  defined  in  Eq.  (2)  hut  hy  using  the 
X y z 

basic  definition  of  angular  momentum  it  is  easy  to  show  that  they  are  equal 
to  the  components  of  the  total  angular  momentum  of  the  molecule  (in  units  h/2Jf), 
including  the  contributions  arising  from  any  internal  rotation  of  the  t<^. 
Similarly  it  can  be  shown  that  p is  the  total  contribution  of  the  motion  of 
the  top  atoms  to  the  z component  of  angular  momentum,  including  both  the 
internal  and  overall  motions. 

This  classical  Hamiltonian  becomes  a qviantum  mechanical  operator  or 
alternately  a matrix  by  regarding  P^,  F^,  and  p as  operators  or  matrices. 

The  commutation  rules^  are 

^ 0.  Klein,  Zeits.  f.  Physik  5§.b  750  (1929). 

PjPj^  - Pj^P,  » -iP^,  (J,k,  £ = x,y,z  in  cyclic  order)  (5) 

and 

pPj  - PjP  = 0 Jo  x,y,z  (6) 

The  first  rule  is  the  stcmdard  one  for  the  conponents  of  angular  mcs&enta  along 
molecule -fixed  axes  while  the  second  follows  from  the  fact  that,  as  an  operator 
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p - - iO/3a)^^_^  (7) 

irheraaB  the  do  not  Involve  a. 

By  inserting  in  E tlse  expressions  for  the  in  terns  of 
and  p^  and  in  ttum  treating  these  latter  aomenta  as  differential  operators, 
the  SehrMlnger  vave  equation  in  the  varlahles  0,^,^,  o could  he  vritten 
down,  hut  this  is  not  necessary  for  the  solution  of  the  prohlesi. 

To  siaplify  the  later  equations,  let 


“ 1 

[e  - D(P^^  + Py^  + /(C-D) 

(8) 

in  which 

D » j (A  + B) 

(9) 

Then 

yf'  - 

b(p/  - Py^)  + P/  - dpP^  + fl>®  + Via) 

(10) 

with 

h = 1 (A  - B)/(C  - D) 

d - 2C/(C  - D) 
f - F/(C  - D) 

(11) 

V - V/(C  - D) 

Then  if  .Alls  an  eigenvalue  of  the  energy  is  given  hy 

W » J(J  + 1)D  + (C  - r)-A-'  (12) 

yf  *f  as  well  as  E,  is  diaigonal  in  J and  M so  only  one  JN  block  need  he 
considered  at  a tlae. 


Synmetry  Considerations 

It  is  well  known  that  the  rigid  asysanetric  rotor  vave  equation  is  liivarlat 
under  the  grovp  of  l80*  rotations  about  each  of  the  principal  axes.  These 
operations  variously  change  the  signs  of  the  P^,  which  leaves  the  energy 
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expression  unaltered.  In  the  j resent  proLlesi  the  or  term  in  pP  not 

* z 

Invariant  unless  p changes  sign,  when  P does.  If  the  potential  energy  V(a) 

z 

is  an  even  function  of  a,  then  the  four  cper-ations  E,  C , C , v'3  whose  effects 
are  listed  below  will  still  lec.'re  H invariant  and  will  form  a group  (the 
”four-gro\^" ) as  before. 


Table  I.  Effect  of  Four -"roup  Symmetry  Operaticrs 


E 


C 


X 


C 

z 


P 

z 


p — >-p 

X;,y  x,y 

I > P 


(a 


a) 


(0 — >0)  (o-^TT-O) 

('f— »<P)  Cf—^Tr+^) 

(%-)■%)  0c-^an-3C) 


■q) 


(«- 


.TT-o’ 


(<r-^Tr+f) 

(pC-^TT-  -)Cj 


I'lp — >4>) 

()6^-n-+  x) 


Consequently  every  solution  of  the  wave  equation  will  belong  to  one  of 
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the  four  species  A,  B ; B , B of  this  group'  as  in  the  case  of  the  eisyamietric 

X y z 

7 

See  G.  W.  King,  R.  M.  Hainer,  and  P.  C,  Cross,  J.  Chem.  Phys. 

27  (19^3). 

rotor,  but  here  the  part  of  the  function  involving  a must  be  included.  Further, 
any  sec\ilar  equation  or  energy  matrix  can  be  factored  into  at  least  four  factors 
by  using  the  above  symmetry. 


For  , in 


Often  V(a)  will  poseesB  addition*.!  syaanetxy  pr^i'perc-ies.  For  exao^'/.e,  in 
nitromethane.  It  presumably  is  invariant  rnder  a — >a  f £'TTk/6j.  where 
k = 0,1>2,3,4,5.  In  addition,  a— >(2TTk/6)  - c;  shovsA  also  le.ave  v un-^hanged. 
These  operations  are  isomorphous  with  t'ae-pcin.t  grrcr^  this  group  is 

not  here  used  in  the  same  way  as  when  gr-yup  theC'?;"  is  ai<plie3.  tc  the  viorat-on 
problem. 

The  symmetry  of  V(o!)  cariX',ot  be  Indis  c-rim-lnfetely  applied  t.c  E because  p 
appears  in  the  cross  term.  Ecwever,  the  "rotaticnje"  a -f  STTk/s  {here 

s =6)  do  leave  H invariant.  So  do  the  ’'reflections”  cx  — ^(2)(k/s)  - a if 
simultaneously  P,  — > -P  . The  four-group  operatlcns  l;!,?ted  earlier  and  the 
8 internal  rotations  a— »G  + 277 k/s  genersd^  ft  <f  % C^erations. 

Table  II  shows  the  resiolt  when  s = 5,  in  SsJfe  is  isomcrphous 

with  the  group 


Table  II.  Character'  Table  for  C-ase  s.  3. 
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The  ons-dinenslonal  species  are  here  labelled  so  that  they  correspond  to  the 
species  for  the  subgroi^  V.  The  degenerate  species  would  become  A 
In  the  subgroup  while  Eg  would  yield  ■¥  B^. 

The  case  s » 6 Is  easily  constructed  from  the  abore  table  because  the 
new  groiqp  Is  obtained  by  Introducing  the  additional  operation  which 
causes  a — + TT.  This  commutes  with  all  operations  above  so  the  new  group 
has  twice  as  many  classes  and  twice  as  many  species  as  the  old.  See  Table  IZX. 


Table  III.  Character  Table  for  Case  s ■ 6. 


9 

These  higher  symmetries  permit  further  factoring  of  the  secular  equation 
or  energy  matrix  if  the  expansion  functions  are  chosen  to  have  symmetries  in 
accord  vith  the  varioiis  species. 

The  dipole  moment^  for  this  model  wooild  ordinarily  be  along  the  z axis . 

If  this  is  the  ceu3e>_^  will  have  the  symmetry  in  Table  II  and  in  the  case 
8 « 6 it  will  belong  to  B . Consequently,  for  this  case,  the  selection  rules 
for  dipole  absorption  will  be  (s  = 3)s  A<— ^B^,  B^^— ^B^.,  — ^®1^  ^2^ — ^ ®2* 

The  rules  will  apply  for  the  case  a b 6 with  the  additional  condition  that  e 
species  go  only  into  e species,  and  o species  into  o species. 

Nuclear  permutation  effects  will  be  governed  by  the  symmetry  also.  The 
subgroip  which  governs  the  exchange  of  the  oxygen  atoms  in  nitromethane  or 
the  fluorines  in  CH^BFg  consists  of  S and  since  to  exchange  these  atoms 
reqxilres  that  a^^ct  + 'Tf ; 7T , which  is  C^Cg.  A^,  B^^,  B^^, 

By^,  Eg^  are  even,  the  other  species  odd  to  this  operation.  Therefore  only 
the  first  set  can  occur  in  nitromethane  (with  ^^0).  In  CB^Fg  these  species 

wi3J.  have  three  times  the  fluorine  spin  weight  of  the  others.  The  operations 

2 

C^,  exchange  two  pairs  of  E atoms  in  the  above  molecules.  All  the  non- 
degenerate species  are  of  species  A in  this  sxibgroip,  all  the  degenerate  simcies 
will  be  degenerate  in  the  subgrov^.  The  two  hinds  of  levels  will  then  have 
equal  proton  spin  weight. 

The  Energy  Matrix 

If  the  asymmetry  (b)  wore  zero  and  there  were  no  barrier  (V)#  the 
reduced  Hamiltonian  3/^would  be 

2 2 
t/  - dpP^  ^ fp 


(15) 
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which  yields  a diaigonal  matrix  with  the  basis  functions 

ilk) 

in  which  S is  the  0,f>  factor  of  the  symmetric  rotor  wave  function  and 

K = 0,+l,+2,  +J  j m = 0,+l,j^2,  ...  oo  . (15) 

These  functions  may  be  used  to  set  vq>  a matrix  for  the  general  form  of 
0^'f  i.e.  the  true  wave  functions  may  be  esqpanded  in  terms  of  the  fiuxctions 
above.  The  esynane  'wCi'ui  in  b will  give  off-diagonal  elements  the  same  as 
those  which  occur  with  the  rigid  asymmetric  rotor  (with  C calculated  for 
framework  only).  The  barrier  potential  V is  usually  assumed  to  be  of  the  form 

V - i-  (1  - cos  8 a)  - I - I (e^®“  + e‘^®“)  (l6) 

in  which  8 is  the  number  of  eqxiivalent  minima  and  (C-D)V^  = is  the  barrier 
height.  The  barrier  will  also  Introduce  off-diagonal  elements.  The  matrix 
for  will  then  have  the  non -vanishing  elements® 

g 

See  ref.  7 for  the  matrix  elements  of  P , P , and  P . 

X'  y'  z 


'^Km,Km.  Km,Kn 


^ - b 

K,m;K+2,m  “ K,K+2 


K,mjK,m+s 


-rv 
k o 


(17) 


Here  the  constant  ^ has  been  Incorporated  so  that  the  eigenvalues  of 
are  related  to  those 


-A.'  + 


(18) 
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Unless  the  harrier  height  is  zero,  this  form  of  ^ corresponds  to  an 
infinite  secular  equation.  Note^  however,  that  there  are  no  elements  connect- 
ing even  with  odd  K values  so  that  the  secular  equation  factors  into  one  for 
even  and  one  for  odd  K values.  (This  is  part  of  the  factoring  into  different 
symmetry  species . } 


Case  of  Free  Internal  Rotation 

When  the  hairier  height  is  zero,  ^ ^ = 0 and  the  energy  matrix 

becomes  diagonal  in  m.  The  secular  equation  therefore  factors  into  one 
block  for  each  value  of  m,  as  veil  as  into  even  and  odd  K factors. 

The  secular  equation  for  a given  veLLue  of  J,M,m  has  the  elements 

(KiK)  « - dmK  - A J ^ 

(KjK  + 2)  = - (K+1)^J[(J+1)^  - (K+l)^jj^  (19) 

The  energy  W is  then  related  to  the  roots  A by  the  eqxxation  (see  Eq,  12), 

W « J(J  + 1)D  + j + Fm^  + (C-D)A  (20) 

9 

This  secular  equation  is  similar  to  tbs  Wang  equation  for  the  ordinary 

^ S.  C.  Wang,  Phys.  Rev.  j4,  24j  (1929). 

rigid  asymmetric  rotor  except  for  the  added  term  -dmK  on  the  diagonal.  This 
term  spoils  the  additional  factoring  possible  in  the  rigid  case. 

As  an  example,  consider  the  case  J » K even.  This  factor  becomes 
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l6  + 4md-X  6/Tb  0 0 0 

6\/5'b  4 + 2 md  -X  ^io  b 0 0 

0 V^IO  b -A  \/^0  b 0 

0 0 V^210  b 4-2  md  - A 6 \fS' \> 

0 0 0 6 yj'b  l6  - 4 md 

For  each  non -vanishing  value  of  ) m | , there  will  be  two  identical  such 
equations,  so  all  levels  are  doubly  degenerate  unless  m = C.  In  the  latter 
case,  the  ordinary  rigid  rotor  levels  are  obtained,  except  that  the  moment  of 
Inertia  about  the  z axis  is  the  moment  of  the  framawork  group  (I^  - I^)  as 
already  seen. 

The  energy  levels  in  the  free  rotation  case  must  of  course  conform  to 
the  symmetry  restrictions  given  for  the  general  case.  Thus,  with  s = J or  6, 
if  m is  not  a multiple  of  the  symmetry  is  one  of  the  degenerate  £ species. 

If  m is  a multiple  of  3,  the  levels  Involve  A or  B species  which  would  be  split 
Into  non-degenerate  components  If  there  were  a sttfficient  barrier.  If  K is 
even,  the  species  may  be  A,  B or  E.;  if  K Is  odd,  B , B or  E„.  If  s ■ 6, 

Z X X ^ c 

the  species  further  divide  Into  even  or  odd  (denoted  by  subscripts  e and  o) 
according  as  m Is  even  or  odd. 

The  selection  rules  for  dipole  reidiation  involve 

dkJ  = 0,  +1  (22) 

as  for  the  rigid  asymmetric  rotor.  The  symmetry  requirements  discussed 
earlier  must  also  be  met.  along  the  z axis,  K cannot  change  parity. 


The  rule  for  m Is 


^ s 0 


(25) 
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This  follows  from  the  fact  that  the  dipole  moment  does  not  depend  on  a 
whereas  the  wavefxmction  Involves  a only  through  the  factor 

Ima 

e 

For  small  asymmetry  the  energy  levels  may  be  eaqpanded  in  powers  of  b. 

For  large  values  of  the  quantum  number  n,  the  diagonal  elements  of  the 
seculeu:  equation  will  be  large  and  the  first  terms  in  the  esqpanslon  in  powers 
of  b may  be  adequate.  Second-order  perturbation  theory  then  gives 

A - - md  K + (b^/8)|j^  - (K-l)^]  [( J+l)®-(K-l)^ /(2K-md-2) 

-(b^/8)|^[j^  -(K+1)^][(J+1)^-(K+1)^ /(2K-md+2)|  + ...  {2k) 

For  very  large  values  of  m, 

X > - md  K.  (25) 

For  transitions  with  s AK  0^  ^ « 0;  the  transitions  with  large  values 
of  m would  then  tend  to  converge  (from  both  sides)  to  a beuid  head.  Of  course 
the  intensities  of  these  higher  members  would  tend  to  decrease  because  of  the 
unfavorable  Boltzmann  factors. 


Case  of  Low  Barrier 

If  the  barrier  height  is  small  but  not  negligible  (conipared  with  the 
quauitlty  F),  a \xseful  solution  can  be  obtained  with  the  Van  Vleck  perturbation 
method. The  reduced  hamiltonian  matrix  of  Eq.  17  is  split  into  am 

See  E.  C.  Kemble,,  The  Fundamental  Principles  of  Quamtum  Mechanics, 
McGraw-Hill,  New  York,  1937^  p.  ’^9^* 


unperturbed,  completely  diagonal  matrix  with  diagonal  e].ements 

Km,Kin 


= - dmK  H 


(26) 
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and  a perturbat 


ion 


with  elements 


K,m;K+2,m"  K,K+2  ? 


K,mjK,iE+B 


- ff ■ - f 


The  Van  Vleck  transformation  reduces  the  elements  off-diagonal  In  m to  order 
5 If  these  are  then  neglected^  the  secxilar  eqtiatlon  becomes 

(K|K)  . - diiK  <•  (l/s)  |(l/jKd-(ani+B)f])  - (l/[Ka-(2m-s)fJ)| 

(K|Kt2)  - 


This  secular  equation  is  similar  In  Its  properties  to  that  for  free  rotation. 
It  is  not  valid,  however,  for  the  case  m = + ^ s,  s even,  because  then  a near 
degeneracy  occurs.  For  this  case,  the  unperturbed  states  with  m * ^ ^ s must 
both  be  taken  Into  consideration.  The  equation  then  has  the  non-vanishing 


elements 

(|.k(|,K)  * - I d S i + (l/s)jTj^/(Kd  - 2sf)J  - A 

(-|,k|-|,K)  - + I d 's  K - (1/s)  jj^/(Kd  + 2sf)^-A 

(|,k|-|,k)  » 5 (29) 


(g->K|g,K+2)  = 

(.|,K|-|,Kt2),= 
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Becavise  of  the  synmetry  (^,K  +g-,K' ) = -K)  this  larger  secular 

eqiiatlon  can  he  factored  into  two  factors^  similar  to  the  way  In  which  the 

Wang  eqiiation  is  factored*  This  leads  to  a splitting  of  the  degeneracy. 

When  m is  not  a multiple  of  g s,  the  levels  are  inherently  douhly  degenerate 

for  all  harriers.  For  m a multiple  of  ^ s the  degeneracy  is  altimately  split 

at  hl]^  enough  harriers.  With  the  approximation  given  here  this  splitting  will 

not  appear  except  for  m = ^ s.  The  higher  the  multiple  of  s the  higher  the 

order  of  perturbation  required  to  demonstrate  the  splitting. 

The  quantum  number  m is  no  longer  a good  quantum  number  when  there  is 

a harrier  hut  if  V is  small  it  can  still  he  used  to  label  the  levels.  The 
o 

selection  rule  ^ = 0 is  no  longer  exact  emd  must  he  replaced  by  the  rigorous 
symmetry  s>?lection  rules  given  earlier. 

High  Barriers 

If  the  harrier  is  sufficiently  high  so  that  the  torsional  levels  of  a 
given  symmetry  are  widely  spaced  compared  with  the  rotational  levels,  a 
different  pertxirhation  treatment  is  possible,  in  which  the  asymmetry  and  the 
coupling  term  are  treated  as  perturbations  by  the  Van  Vleck  procedure.  This 
leads  as  before  to  a rotational  secular  equation.  The  unperturbed  operator 
for  the  energy  is 

H = DP^  + (C-D)  + Fp^  + V(a)  (30) 

o z 

while  the  perturbation  ojjerator  is 


= I (A-B)(P^  - Py)  - 2CpP^ 


(31) 
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can  be  diagonalized  vlth  the  basis  functions 


(32) 


In  vhlch  symmetric  rotor  functions  are  used  as  before  and  the  U's  are  eigen- 
functions of 


+ V(a)l  U . 


, , ^ 22  E \*TT  K A 

Li  'vK 
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The  torsional  states  are  described  by  the  quantum  numbers  v)^,  where  v Is  the 
principal  quantum  number  of  the  vibrational  levels  Emd  Is  a degeneracy  Index. 

For  the  potential  function  of  equation  (16)^,  the  non -degenerate  = O) 
eigenvalues  of  the  torsional,  equation  (j>3)>  which  Is  now  related  to  the 
Mathleu  equation,  can  be  obtained  fr<xn  publlshed^^  tabulations  through  the 


Tables  Relating  to  Mathleu  Functions,  Columbia  University  Press, 
Hew  York,  1951. 


coxmectlon 

®vo  = f ^ I ^ FbOy(S)  (54) 

in  which 

S - (4/s^)(vyF)  (35) 

and  the  elgenvedues  be^  and  bo^  are  given  in  the  tables  as  functions  of  the 
parameter  S.  The  values  of  the  quantum  number  v = 0,1, 2, 3, 4,  etc.  are 
Identified  respectively  with  the  eigenvalues  be^,  bOg , beg,  boj^^, 


IT 


The  degenerate  levels  have  not  been  tabiilated  but  can  be  calculated  to  any 
desired  degree  of  accuracy  by  expansion  of  the  eigenfunction  in  a Fourier  series. 
Another  method  is,  however,  easier  and  is  reasonably  accurate  for  high  barriers. 
Make  the  substitution 


Us.  - ufa) 

ViV  ' ’ 


in  Eq.  (33).  where  if  = 0,  1,  or  -1.  The  equation  for  u is  then 

+ 2Iifcp  + fJ?^  + vJ  u = E^^u  (58) 

If  if  = 0,  the  eqviatlon  gives  the  non -degenerate  elgenvedues  smd  their  eigen- 
functions, both  of  which  are  available  as  tabulations.  If  if  = ^1,  the 
degenerate  levels  are  obtained.  Bnt  the  terms  involving  in  Eq,  (58)  can 
be  considered  as  perturbations  on  the  if  = 0 case,  if  the  barrier  is  hig^h 
enough.  The  term  in  obviously  jiist  subtracts  from  E^.  To  the  second- 
order  the  effect  of  the  terms  in  if  suad  is 


E 


VO 


+ Fif^' + 


/(E. 


VO 


E , ) 

v'o 


(39) 


in  which  the  sum  is  over  all  the  non -degenerate  levels  and 


. r.: 


Ou,.y3o)ao 

r pc 


r pc  ^ "N  , 

21  mA  A , /£  A^  2.  A^ 

vm  v'ny  / vm  ^ v'm 

m=-oo  7 m--oo 


Sere^  we  have  ejqpanded  the  non -degenerate  eigenfunctions  as 

lm!?C 


U (a) 
vo'  * 


-fL 


m=-cc 


A e 
vm 


(40) 
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These  Fourier  coefficients  may  he  found  from  the  tah\ilatlons  of  reference  11 
hy  the  Identification 

V-sk  = Vk  - k^O 

A = 2De^(S)  for  even  v 

T,o  o'  ' 


\,-8k  “ “\,8k  “ ^ 

The  perturbation  of  Eq.  (31)  now  needs  to  be  Introduced  by  the  Van 
Vleck  procedure.  However,  for  the  degenerate  torsional  states,  it  is  con- 
venient to  combine  this  perturbation  with  that  given  above,  yielding  the 
secvQar  equation 

(K|K)  = - d>)te  ^j\r  ^ 

(k|k  + 2)  » ^b‘|'[j^  - (K  + if]  [(J  + 1)^  - (K  + 

with  W » J(J+1)D  + E^^  + + JSl'  (C’-D) 

C’  =.  C(1  + kC^) 

d’  » 2C(1  + 4F^)/(C'-D)  (^^2) 

F’  » F(1  + 4F^) 

/O  « 21^  IP  , |V(E  - E , ) 
f , rw’  ' ' VO  v’o' 

V 

b*  - |<A-B)/(C’-D) 


Very  Hlf^  Barriers 

If  the  barrier  V is  sufficiently  high,  the  lower  energy  levels  should 

approximate  those  of  a rigid  rotor  with  mooents  of  inertia  I , I , I 

X y z 

(not  Ijg  “ ^ harmonic  torsional  vibrator  with  a reduced  moment  of 


inertia 


and  potential  energy 


I (dVaa^)oa^ 

This  result  can  be  shown  to  follow  from  Sq.  (4l)  and  (^2),  since 

4^  > -1  (45) 

for  high  barriers  In  which  harmonic  oscillator  functions  cem  be  used  for 
In  Eq.  (40). 

The  energy  levels  of  rotation  are  then  cej-culable  by  the  standard  methods 
\ised  for  the  rigid  asymmetric  rotor.  The  torsional  levels  will  be  degenerate^ 
three-fold  If  n a 3,  six-fold  If  n » 6.  lowever,  some  of  this  degeneracy  may 
be  eliminated  when  permutation  symmetry  occurs,  as  In  nltronethane,  whose 
levels  will  be  only  three -fold  degenerate. 

Case  of  Small  Asymmetry 

In  cases  for  which  the  approximations  so  far  described  core  all  Inadequate, 

12 

anothar  approach  ceui  be  used,  as  long  sus  the  asymmetry  of  the  molecule  Is 


This  Is  very  close  to  that  Introduced  by  Dennison,  ref.  1. 

small.  Let  the  reduced  energy  matrix  3^'  of  Bq.  (10)  be  set  up  in  terms  of 
the  basis  functions 

In  which  S Is  a symmetric  rotor  function  and  Is  a periodic  function 

of  a which  Is  an  eigenfunction  of  the  equation 

[fp^  - dl^  + V(a)jy^j^(a)  - W^jj‘>^(o) 


(47) 
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with  eigenvalue  K enters  as  a parameter.  The  non-vanishing  matrix 

elements  for  .^^then  become 


^^,yK 

• '5/'  = b A 

v,KiV,K+2  °K,K+2  vK,v'K+2 

(48) 

otr 

where 

\,Kjv',K+2  ” f y'vK^yy' 

(1^9) 

which  enters  because  for  different  K's  are  not  orthogonal. 


The  existence  of  the  terms  off-diagonEuI  In  v spoils  the  factoring  Into 
V -blocks  but  these  terms  are  normally  small  and  can  be  reduced  to  lover  order 
by  a Van  Vleck  transformation.  This  yields  a rotational  secular  eq[uatlon 

(k|k)  - - -A-' 

* l'’K,K+2i  hvK,v'K+2|  " '^v'K+2^ 

* K.,K-2l  ^ t^^v'K=2l 

(k1k+2)  - \K,y/K±2 

'^vK  " ’^v'Kj^ 

15 

Koehler  and  Dennison  have  shovn  that  the  quantities  can  be  obtained 
J.  S.  Koehler  and  D.  M.  Dennison,  Phys.  Rev.  57#  1006  (1940). 
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frost  a continued  fraction  of  the  form 


L_  1 

5 "c  ... 

Ill 

- v5-^-®r5-  ••• 


in  which 


- g)  > 9.  - 0,+l,+2 

g - [(CK/f)  -)^)]/8,'ht=  0,  1,  2,  8-1  (52) 

^nK  “ I ^o  ■ CcV/f(C-D)]  + s^f  ^ 

To  calculate  the  A'e  continued  fractions  may  also  he  used*  The  function 


’VV  is  ejqjanded  in  a Fourier  series; 

**  o<J 


Ti=..0O 


cK,m 


This  is  substituted  in  the  differential  eq'jation  (47)  and  leads  to  a set  of 

simultaneous  equations  for  the  coefficients  A: 

A ^ + G A „ + A « ■>  0 

nK^ui-s  m nK^m  nK^m-fs 


These  simultaneous  equations  yield  the  contintied  fraction 


^nK,iM-8  _ -_1  1 1 

■^hK,m  ~ ®m+s  ” ®m+2s 


The  normalization  condition 


m=-co 


nK,m 
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envies  the  determination  cr  the  A' 3 to  he  cociplated. 

There  vlll  of  course  he  an  infinite  nuaiber  of  roots  for  each  valvie  of  K 
and  VC . These  roots  are  moat  conveniently  calculated  if  none  of  the  denominators 
“3"  ••*  small.  If  »ne  of  these  should  turn  out  to  he  small,  the  continued 
fraction  should  he  transfonaed  so  that  the  offending  appears  in  the  leading 
positionj  i.e. 


•?  - M 1 1 


(55) 


For  each  root  sought,  there  will  be  an  appropriate  choice  of  the  leading  term 

which  will  cause  the  most  rapid  convergence. 

The  varlotis  values  of  ^ yield  the  different  component  levels  which  in 

the  limit  of  high  harriers  eone  together  to  form  a (degenerate)  harmonic 

' t 

oscillator  level.  However,  Koehler  and  Dennison  ^ showed  that  the  roots 
Bxe  periodic  functions  of  K (as  may  also  he  seen  from  S5.  51)  with  period 
Fs/C.  Furthenaore,  if  this  periodic  function  is  known  for  the  case  0, 
it  can  he  used  to  calculate  the  roots  for  the  other  vsilues  of  by  a l^hase 
shift,”  i.e.  hy  reading  off  the  veilue  of  ^ for 

K'  = K - F>f/C 

from  the  same  periodic  curve  used  for  0. 

For  sufficiently  high  barriers  the  dependence  of  ^ on  K and  ^ vanishes. 
This  suggests  eigpandlng  £ as  a Fourier  series  (only  cosine  terms  will  appear) 
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5 = 5o  + 5 ^ ^ 8 


cos 


hTTg 


For  high  enough  harriers  only  a few  terms  are  required.  The  coefficients 

fo«nd  hy  solving  the  continxied  fraction  for  the  necessary 
nuniber  of  veJ-ues  of  g,  conveniently  chosen  to  make  it  easy  to  solve  for  the 
required  number  of  ^'s.  The  value  of  can  also  he  obtained  from  available 
tables  ^ hy  means  of  the  relations 


with  S a kV/e^f  = 4V  /s^ 

o'  o'. 

Then  the  A's  may  he  written 


oo 


K'.nK 


'■ttzL 

jQHa:*00 


'n'K’ 


,m  '^nK,m 


In  sane  eases,  especially  when  h is  8\ifficlent].y  low,  the  (k|k+4)  elements  of 

the  secular  eqsiatlcn  (^O)  can  he  neglected  altogether.  However,  this  procedure 

is  not  always  vsJJ.d.  Under  such  circumstances  the  secular  equation  ceua  he 

made  easier  to  solve  hy  varlovis  transformations  which  reduce  the  order  of  the 

troublesome  K,Kt4  el«nents  so  thnt  they  can  he  ignored  and  the  remaining  equation 

are 

treated  hy  the  continued  fraction  method.  These  transformatlons/slmilar  to 
Van  Vleeh  transfozmatlons  hut  have  to  he  adapted  to  the  several  cases  which 
arise  because  of  near -degeneracies  and  will  therefore  not  he  detailed  here. 


